Lattice model for kinetics and grain size distribution in crystallization 
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We propose a simple, versatile and fast computational model to understand the deviations from 
the well-known Kolmogorov-Johnson-Mehl-Avrami kinetic theory found in metal recrystallization 
and amorphous semiconductor crystallization. Our model describes in detail the kinetics of the 
transformation and the grain size distribution of the product material, and is in good agreement with 
the available experimental data. Other morphological and kinetic features amenable of experimental 
observation are outlined, suggesting new directions for further validation of the model. 



I. INTRODUCTION 

Mechanical, electronic or magnetic properties of many 
polycrystallline materials depend not only on their chem- 
ical composition, but also on the kinetic path of these ma- 
terials toward the non-equilibrium state. Recently, the 
interest on thin film transistors made of polycrystalline 
Si and Si-Ge grown by low-pressure chemical vapor depo- 
sition has been driven by the technological development 
of active matrix addressed flat-panel displays jlj] and thin 
film solar cells, [0 With these and similar applications in 
mind, the capability to engineer the size and geometry 
of grains becomes crucial to design materials with the 
required properties. 

In general, crystallization of most materials takes place 
by a nucleation and growth mechanism: |j| Nucleation 
starts with the appearance of small atom clusters (em- 
bryos). At a certain fixed temperature, embryos with 
sizes greater than a critical one become growing nuclei; 
otherwise, they shrink and eventually vanish. Such a 
critical radius arises from the competition between the 
surface tension, 7, and the difference in free energy be- 
tween the amorphous and crystalline phases, Ag, that 
favors the increasing of grain volume, yielding an energy 
barrier that has to be overcome to build up a critical nu- 
cleus. For a circular grain of radius r, the free energy 
takes the simple form 



AG = 2irrj — nr Ag. 



(1) 



The free energy AG has a maximum, the energy bar- 
rier, at the critical radius r* — j/Ag. Subsequently, 
surviving nuclei (r > r* ) grow by incorporation of neigh- 
boring atoms, yielding a moving boundary with temper- 
ature dependent velocity that gradually covers the un- 
transformed phase. Growing grains impinge upon each 
other, forming a grain boundary, and growth ceases per- 
pendicularly to that boundary. Therefore, the structure 
consists of vertices connected by edges (grain boundaries) 
which surround the grains. The number of edges joined 
to a given vertex is 3. In some cases, at high temper- 
atures these boundaries move until they reach a more 
favorable equilibrium configuration (in two dimensions, 
the equilibrium angles at a vertex are 120°). |J] 



In the past few years, the belief that this picture is far 
too simple to properly describe nucleation-driven crystal- 
lization has progressively spread among the researchers in 
the field. This is chiefly due to two problems: On the one 
hand, this theory of nucleation and growth predicts an 
energy barrier much larger than the experimental one, 
implying that nucleation would be hardly probable at 
available annealing temperatures. Q On the other hand, 
it is known that in crystallization of Si over Si02 sub- 
strates, nucleation develops in the Si/Si02 interface due 
to inhomogeneities or impurities that catalyze the trans- 
formation. H Therefore, a theory of homogeneous nucle- 
ation and growth is not entirely applicable to the referred 
experiments as well as to other examples reported in the 
literature. 

In addition to the difficulties above, it is clear that 
the transformation kinetics is also problematic. It is gen- 
erally accepted that the fraction of transformed mate- 
rial during crystallization, X(t), obeys the Kolmogorov- 
Johnson-Mehl-Avrami (KJMA) model, |s) according to 
which 



X(t) = 1 - exp(-at m ), 



(2) 



where a is a nucleation- and growth-rate dependent con- 
stant and to is an exponent characteristic of the exper- 
imental conditions. Two well-defined limits have been 
extensively discussed in the literature: When all the nu- 
clei are present and begin to grow at the beginning of 
the transformation, the KJMA exponent, to, is equal to 2 
(in quasi- two-dimensional growth like thin films), and the 
nucleation condition is termed site saturation. The prod- 
uct microstructure is tessellated by the so-called Voronoi 
polygons (or Wigner-Seitz cells). On the contrary, when 
new nuclei appear at every step of the transformation, 
to = 3 and the process is named continuous or homo- 
geneous nucleation. Plots of log[— log(l — X)] against 
log(f) (called KJMA plots) should then be straight lines 
of slope to. Although in some cases, the KJMA theory 
explains correctly the transformation kinetics, its general 
validity has been questioned in the last few years, [p|jic|l 
and several papers have been devoted to understand this 
question in different ways, [[ll] |l3| However, there are still 
some open questions: An exponent between 2 and 3 is 
experimentally obtained (between 3 and 4 in three di- 
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mensions); |Q the KJMA plots from experimental data 
do not fit to a straight line in some cases; and, finally, 
the connection between geometrical properties (grain size 
distributions) and the KJMA exponent is not clear. 

In this paper we report on a detailed investigation of 
a probabilistic lattice model which relates in a clear-cut 
way the mentioned problems to the inhomogeneities in 
the sample, i.e., the fact that heterogeneous nucleation 
takes place. Indeed, heterogeneous nucleation is rather 
common in nature due to impurities or substrate cavi- 
ties resulting from roughness, among others. Within our 
model, the connection between such heterogeneous nu- 
cleation and the deviations from the simplest nucleation 
picture become evident. Furthermore, as we will see be- 
low, our model predicts measurable quantities, such as 
the grain size distribution or the KJMA exponent, which 
are in good agreement with the experiments. Our pa- 
per is organized according to the following scheme: in 
Sec. II we introduce our model and discuss in depth the 
relationship between its defining parameters and physi- 
cal ones. Section III collects the results of an extensive 
simulation program which establishes the main features 
of the model. Finally, Sec. IV discusses the connection 
between our model and experiments, and concludes the 
paper by summarizing our main findings and collecting 
some prospects and open questions. 



II. THE MODEL 



Simulation proceeds in discrete time steps of duration 
t. The system evolves by parallel updating according to 
the following rules and considering that initially all the 
material is untransformed, i.e., g(x, 0) = for all lattice 
sites: 

• An already transformed site remains at the same 
state forever. 

• An untransformed potential site may become a new 
non-existing state (i.e., crystallizes) with probabil- 
ity n (nucleation probability) if and only if there 
are no transformed nearest neighbors around it. 

• An untransformed site (including potential sites) 
transforms into an already existing transformed 
state with probability g (growth probability) if and 
only if there is at least one transformed site of that 
type on its neighborhood. The new state is ran- 
domly chosen among the neighboring grain states, 
if there are more than one. 

Note that we have termed g as growth probability and 
not growth rate. The actual growth rate is a non-trivial 
function f(g), because when g < 1 the grains grow with 
a rough boundary. For the model parameters, we expect 
a functional form n ~ e -E„/k B T anc j ^ e -E g /k B T ^ 

where E n and E g are the energy barriers for nucleation 
and for growth, respectively (see below). Hence, temper- 
ature is implicit in the model parameters. We discuss 
these relationships in depth in the next subsections. 



A. Evolution rules 

Our model is based in some previous ideas by Calm 
Jl4],[l5| and Beck, |l6| and its key proposal is that the 
material is not perfectly homogeneous but, on the con- 
trary, it contains regions with some extra energy (regions 
with some order produced during deposition or substrate 
impurities) at which nucleation is more probable. Our 
aim in this section is to provide a detailed description 
of our model (largely expandin g th e preliminary, short 
report presented earlier in Ref. []17|] ) , and how the basic 
idea mentioned above is implemented in it. 

The model is defined on a two-dimensional lattice 
(square and triangular lattices were employed with essen- 
tially similar results) with periodic boundary conditions; 
generalizations can straightforwardly be done to any spa- 
tial dimension. In the beginning (t — 0), every lattice site 
(or node) x belongs to a certain grain or state. We rep- 
resent the situation at x by q(x, t) — 0, 1, 2,. . . , the state 
being that of an untransformed region. The lattice 
spacing is therefore the experimental resolution, usually 
greater that r* . Following the idea that the amorphous 
phase has random regions at which nucleation is favored, 
we choose a fraction c of the total lattice sites and label 
those as able to nucleate. We term these energetically 
favorable sites potential nuclei. 



B. Physically relevant magnitudes 

As we mentioned above, the crystalline fraction is ap- 
proximately given by Eq. (||) , with some exponent m de- 
pending on the dimensionality and type of nucleation. 
Experimentally, the crystalline fraction is measured from 
the intensity of the peaks of X-ray diffraction of the mi- 
crostructure as material transforms from the amorphous 
to the polycrystalline phase. In the following, we will 
assume that there is not any preferential direction, |lg| ] 
that is, n is the same for all potential sites, and g is the 
same for all grains. 

The other experimentally measurable magnitude is the 
grain size distribution, P(A), defined as the fraction of 
grains with a given area A. To compare with simulation 
results, we will usually plot the normalized distribution 
of reduced area A' = A/ A, where A is the mean area: 

/>OC 

A= / AP(A)dA. (3) 
Jo 

This distribution changes dramatically with nucleation 
conditions. JT~g| ] Some of the available experimental data 
are given in terms of the distribution of grain diame- 
ters, P{d). As we will demonstrate below, this distribu- 
tion is equivalent to the distribution of effective diameter 
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A 1 / 2 /it (or simply A 1 / 2 ), which is computationally less 
expensive to calculate. Hence, we will present our results 
in terms of the effective diameter. 



C. Time and length scales 

To begin with, let us show that the potential sites, dis- 
tributed randomly throughout the system, define a char- 
acteristic length given by the probability distribution of 
nearest neighbors. Suppose we have N randomly poten- 
tial sites in a L x L system. The mean concentration 
of potential sites is c = N/L 2 . We may ask about the 
probability of finding a number k < N of these sites in a 
region of area A. This probability is given by the bino- 
mial distribution: 



Pk(N,A) 



p k (l-p) 



N-k 



(4) 



where p = A/L 2 . 

Taking the limit L — > oo, N — > oo, while keeping 
N/ L 2 = c, the expression (J4j) tends to the Poisson prob- 
ability distribution, 



P k (A) = i^fe-^ 



(5) 



If, as stated above, we suppose that the system is 
isotropic, we may write the probability of finding k po- 
tential sites in a circle of radius r as: 



2„\k 



Pk{r) 



(ixr c) 



(6) 



So, if Eq. (||) is the probability of finding n potential sites 
in a disc of radius r, then the probability of finding no 
grains is 



Pair) 



(7) 



and the probability of finding at least one neighbor at a 
distance less than r is 1 — Po( r )- This is precisely the 
probability distribution of nearest neighbors. In other 
words, we can obtain the probability density of finding 
at least one neighbor between r and r + dr as follows 

pMdr = —(l- PJr)) dr = 2tt rce'^ dr. (8) 
dr V / 

The first moment of the distribution is the mean distance 
among potential sites 



d m = I rp(r)dr = c 



= .-1/2 



(9) 



On the other hand, the grains grow with constant ve- 
locity. For defmiteness, let us take the growth probability 
g to be 1; we will see below that the results of simulations 
for other values of g can be reproduced from simulations 



with g = 1 conveniently rescaled. With this choice, the 
grain radius grows according to the law r(t) = fit, where 
fl is a geometrical coefficient that depends on the un- 
derlying lattice. Thus, we may define the mean time at 
which the growing grains will impinge, or overlap time, 
as flt = c -1 / 2 or, in general, i.e., ignoring the details of 
the lattice, t ~ c -1 / 2 . 

The characteristic time scale arising from the concen- 
tration of nucleation sites is not the only one: Indeed, 
the nucleation probability defines another characteristic 
time. Being more specific, the number of sites that have 
nucleated per unit time is proportional to the available 
ones 



dN(t) 
dt 



= n[N max - N(t)}, 



where N„ 



cL . Thus, we have 



N(t) = 7V max (l - «.—*) => p{t) = c (l - e- nt ), (10) 

p{t) being the concentration of already nucleated poten- 
tial sites at time t. In view of this, we define the charac- 
teristic nucleation time t n — l/n. As we will see below, 
the competition between time scales characterizes the fi- 
nal microstructure. 

In a general case, some of the potential sites will be 
covered by other growing grains and therefore their nu- 
cleation is inhibited. The mean distance of the potential 
sites that become actual grains is, replacing c by p(t) in 



d m (t) 



1 



-1/3(1 - 



(11) 



If t n <C to, almost every potential site nucleates before 
grains impinge upon each other. We term this situation 
fast nucleation, and in terms of our model parameters 
it means that n 3> c 1 / 2 . This situation is similar to site 
saturation nucleation, in which every potential site nucle- 
ates at t = 0. The KJMA exponent will be close to 2 and 
the grain size distribution will be similar to that of site 
saturation. Note that, when n = 1, the exact limit is ob- 
tained for every concentration c < 1, but concentrations 
c close to 1 yield a mean grain size of just a few times 
the critical radius, r* , which in fact has not much to do 
with the experimentally measured values. In this case, 
t n is approximately equal to the simulation time step, r, 
so the characteristic time scale is Tf ast ~ t a ~ c -1 / 2 . 

Analogously, if t n ^> t D then c 1 / 2 3> n and growing- 
grains will overlap potential sites before these have nucle- 
ated, forcing the number of nucleating grains to decrease 
with time. As new grains still appear at every stage of the 
transformation, we expect approximately homogeneous 
nucleation, and correspondingly a KJMA exponent close 
to 3. We term this situation slow nucleation. Compar- 
ing the radii of the grains with the mean distance among 
them we find the characteristic time of the process: 



3 



1 ^ 1 

c x /2[l - ex p(-nr slow )] ~~ (cnr slow ) 1 /2 

~ ^(^slow) — ^^"slowi 



and hence 



''"slow 



(cn) 1 / 3 ' 



(12) 



The important point, however, is the fact that between 
both limits we will find a wide range of KJMA expo- 
nents and grain size distributions, consistently with the 
experimental results. 



III. NUMERICAL RESULTS AND DISCUSSION 



A. Isolated grain shapes 



An isolated grain, i.e., a grain completely surrounded 
by untransformed material, grows isotropically. Thus, in 
a continuum medium, the grain boundary is nearly a cir- 
cumference. Nevertheless, the shape of such propagating 
interfaces in our model depends strongly on the underly- 
ing lattice. For example, in the limit case g = 1, a grain 
growing in a square lattice is square shaped, whereas if 
growing in a triangular lattice it is hexagonal shaped. As 
the growth probability g diminishes, the underlying lat- 
tice effects seem to vanish, and grains are approximately 
circular, with a rough boundary. In Figs. |l| and |^ wc 
show the dependence of the grain shape on the growth 
probability, varying g from 0.1 to 1, on square and tri- 
angular lattices respectively. We see that for g < 0.4 
the shape of an isolated growing grain becomes practi- 
cally independent of the lattice, whereas for larger values 
of g, the grain shape exhibits the influence of the lat- 
tice geometry. It is important to note that this does not 
occur when many grains grow simultaneously, as in this 
case the grain geometry is determined by the succesive 
impingement with its neighbors. 

In connection with the last remark, it is interesting to 
consider another issue related to boundaries, namely that 
of boundaries between different grains. Let r% and r 2 be 
the radii of two circular grains; their boundary is then 
defined by the equation JL9| | 



B. Kinetics 

We have simulated 1000 x 1000 triangular and square 
lattices and averaged the outcome of 50 different realiza- 
tions for each choice of parameters (characteristic simu- 
lation times are about 15 to 45 minutes in a Pentium II 
personal computer). The crystalline fraction ranges from 
to 1, so wc define the typical simulation (or experimen- 
tal) time as the time t 1 / 2 at which X(t 1 / 2 ) = 1/2. As a 
check on our ideas, we begun by verifying the dependence 
of this parameter on the time scales defined above. In 
Figs. |] and g, we plot t 1 / 2 for different parameters in the 
fast and slow nucleation limits. A very good agreement 
is observed with the expected beha vior of 1 i/ 2 ~ Tfast and 



1/2 



T s iow discussed in Sec. II C. Therefore, we can be 



confident that the expectations drawn above about the 
behavior of the model, based on theoretical considera- 
tions, will be fulfilled. 

The first key feature to analyze relates to the crystal- 
lization kinetics as seen through KJMA plots. Our re- 
sults show that those are not the straight lines predicted 
by the KJMA model: This can be best seen by looking 
at the transient KJMA exponent, defined as 



m{t) 



d(logi) 



pog(-log(l -*))], 



(14) 



ri +V g ti = r 2 +v„t 2 , 



(13) 



Figure || shows that the KJMA exponent always de- 
creases from its initial value to an asymptotic, time inde- 
pendent one; correspondingly and in agreement with the 
experiments, KJMA plots approach straight lines only at 
late times. We note that, in determining m(t), care has 
to be taken from the computational point of view as in 
some cases the number of steps needed to complete the 
transformation is too short. In addition, it is necessary to 
remove the last few instants of the time evolution, as they 
exhibit large finite size effects. The asymptotic value is 
the one we take from simulations and the one plotted in 
Fig. [7] showing the dependence of the KJMA exponent 
with the potential site concentration c. Alternatively, 
Fig. |8| depicts the dependence of KJMA exponent on the 
nucleation probability n. We thus see that there is a 
large variability of the KJMA exponent, covering all the 
range between 2 and 3 in this two-dimensional case, that 
depends on the relationship between the nucleation prob- 
ability n (i.e., the nucleation rate) and the concentration 
of nucleation sites c. This result is a step beyond KJMA 
theory, and agrees with the fact that experiments offer 
very different results, with exponents between 2 and 3. 



where v g is the growth velocity and t\$. the elapsed time 
since each grains nucleated. When grains started to grow 
at the same time {t\ = t 2 ), the boundary is a straight 
line. Otherwise, it is a hyperbola. In Fig. || we plot 
two examples of interfaces in which, in spite of the fact 
that interfaces are noisy, both characteristic curves are 
revealed. 



C. Grain area and grain diameter 

In order to further check the model results, we have 
to compare the grain size distributions with some well- 
accepted theoretical ones. Although these distributions 
are obtained phenomenologically, the agreement with ex- 
periments and simulations is very good. Under some 
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assumptions about the mean number of neighbors of a 
nucleation center, Weire et al. proposed a simple distri- 
bution for site saturation [E0| 



P(A>) = {A') 



'exp[-aA']/T(a), 



(15) 



where a ~ 3.65, and A' = A/ A is the reduced area. In 
Fig. ^ we plot the normalized grain size distribution (cir- 
cles) for different parameters for which to ~ 2, i.e., site 
saturation, and compare it with Eq. ( |l5| ) (solid line). 

Similarly, in the case of homogeneous nucleation, a 
simple (but not so accurate) expression has been pro- 
posed §IJ 



P(A') = exp[-A']. 



(16) 



Our model shows some slight deviations from this equa- 
tion, as seen in Fig. |[ Interestingly, these are the same as 
in other model simulations, |2lf| and, in addition, we have 



to keep in mind the applicability limitations of Eq. (16). 
plf Therefore, we believe that the behavior displayed by 
our model is also fully satisfactory in this limit. 

Once we have checked the validity of the model in the 
well-known limits, we report on the influence of the nu- 
cleation probability, n, and the potential site concentra- 
tion, c, on the grain size distribution. In Fig. [l0| we 
plot several grain size distributions when we pick both 
parameters along a line going from the slow to the fast 
nucleation limit. In so doing, we cross from a extended 
distribution to a stretched one, as we would expect in 
view of Eqs. (|l|) and @. 

Let us now turn to the issue of the mean grain size. 
As we have pointed out, in the fast nucleation limit the 
characteristic length scale is related to the mean poten- 
tial site distance, c -1 / 2 : In this case, we expect the mean 
grain diameter to be proportional to that scale. In Fig. O 
we show this linear dependence of the mean area on 
On the contrary, in the slow nucleation limit, when the 
concentration c is relatively large, the grains grow on an 
effective homogeneous medium. Roughly speaking, the 
mean distance among potential sites is so small that the 
grain radii is very soon larger than this distance. Thus, 
the characteristic length scale, d, is that of the grains 
when they impinge upon each other. As the grain ra- 
dius grows linearly with time, we expect d ~ tj/2, so 
A 1 / 2 ~ (crt) -1 / 3 and A ~ (cn)~ 2/3 . Figure [jj] confirms 
that this simple analysis is very accurate. 

Finally, there are two questions we announced in Sec. 
II whose validation has been left postponed. We now ad- 
dress these points, beginning by that of the effect of the 
parameter <?, which so far we have restricted to g = 1. 
For every value of g, the growth rate, f(g), defines a char- 
acteristic time related to the temporal scale at which the 
grains spread on the amorphous substrate. Thus, we ex- 
pect that by rescaling the simulation time step r — » f(g)r 
(with f(g) — » f2, as g — » 1, Q, being the geometrical co- 
efficient introduced in Sec. IIC), the mean grain size will 
depend only on the ratio n/f(g). We have not been able 



to obtain an analytical expression for f{g) but we can 
calculate it numerically for the required <?, by growing an 
isolated grain. In Fig. [l^ we show the excellent collapse of 
different effective diameter distributions for several cou- 
ples (n,g) with constant n/f(g). This result shows that 
the outcome of the simulations reported here for g = 1 
truly represents, except for a factor, the model charac- 
teristics for other values of g. 

The other pending question is related to the mean 
grain diameter. So far, we have discussed our results 
in terms of the mean grain area or the mean effective 
diameter size. To verify whether the effective diameter 
distribution is the same as the real diameter distribu- 
tion, which is computationally much more demanding, 
we have compared them in several cases. The compar- 
ison is shown in Fig. [l4|, by plotting the referred nor- 
malized distributions. The correlation between both sets 
of points is greater than 99.9%, allowing us to conclude 
that the reports above in terms of areas carries over to 
the mean diameter picture without significant changes. 



D. Mean number of neighbors 

Some theoretical approaches to equilibrium crystal- 
lized configurations deal with the mean number of neigh- 
bors, N nn , |2l[] or equivalently, considering the final prod- 
uct as a polygon tessellation of space, the mean number of 
sides of those polygons. If the material is divided in equal 
size hexagons, this distribution is P(N nn ) = 5(N nn — 6). 
In Fig. [l5] we plot the numerical distribution of near- 
est neighbors for site saturation and homogeneous nucle- 
ation. The asymmetry and the variance of the mean 
number of neighbors are the main differences in both 
limits. The inset in Fig. |l5| shows the mean number of 
neighbors and the corresponding changes in variance for 
different parameters. Clearly, the distribution spreads 
out and loses its simmetry in homogeneous nucleation. 
Furthermore, computing the mean number of neighbors 
against the nucleation time for all of the grains in the 
sample we find that the younger grains have less number 
of sides than the older ones, which explains this asym- 
metry. Hence, this distribution can be another element 
of comparison with experiments. We remark that sec- 
ondary crystallization (or abnormal grain growth) is due 
to these deviations from the ideal configuration. 



E. Temperature and applicability of the model 

To conclude our analysis of heterogenous nucleation, 
we present some results of the influence of temperature 
in product properties. In addition, this will allow us to 
show that the model gives consistent results when real- 
istic parameters are chosen to reproduce an actual ma- 
terial. The mean grain area in homogeneous nucleation 
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of two-dimensional disks is given by the simple relation: 



A 



Go 



2/3 



(17) 



where No and Go are the nucleation and growth rates 
respectively. Identifying iVo with n and Gq w ith /(g), we 
have obtained similar results (see Sec. IIIC). As nucle- 
ation and growth are activated processes, we postulate 
an Arrhenius-like dependence of nucleation and growth 
probabilities: 

n~exp[-E n /k B T\, f(g)~exp[-E g /k B T}. (18) 

In homogeneous nucleation, as we have reported, we can 
redefine n and g to set g = 1; hence, the temperature 
is introduced in our model by means of the nucleation 
probability 

n -> ri = n/f{g) = n' cxp[-(E n - E g )/k B T], (19) 
and 5 = 1. 

As an example, if we want to model nondendritic 
Si crystallization, we may use the experimental activa- 



tion energies |22j]: E n = 5.1 eV and E g = 3.2 eV. 
Then, A ~ exp[E a /k B T], where from Eq. @) E a = 
2{E n - E g )/3 ~ 1.27 eV. In Fig. |6| we plot the mean 
grain size vs. 1000/T. The slope gives E a = 1.26 ± 0.01 
eV, which is consistent with the introduced values. Thus, 
the model provides a simple tool to analyze crystalliza- 
tion experiments: Setting the activation energies as the 
program input, we just have to choose a realistic value of 
n Q (e.g., in terms of the final number of grains) and tune 
the degree of heterogeneities, c, in order to compare with 
the experiments. 



IV. CONCLUSIONS 

As we have seen, the model proposed in this paper 
provides very accurate and detailed spatial and temporal 
information about the system evolution: Crystalline frac- 
tion, mean grain area, KJMA exponent, or mean number 
of neighbors. The main features observed in experiments, 
such as non-integer KJMA exponents or different types 
of grain size distributions are very well reproduced by 
the model. We must conclude, then, that the model cap- 
tures all the physical ingredients involved in the crystal- 
lization process: In particular, it points out to the in- 
homogeneity of the nucleation phenomenon (which can 
arise because of the structure of the amorphous mate- 
rial itself, or because of defects at the substrate-material 
interface, for instance) as the key feature governing the 
crystallization kinetics and the resulting grain textures. 
In view of this, we propose this model, very unexpensive 
in terms of computing time, as a versatile way to incor- 
porate other physical ingredients as boundary migration, 



preferential grain growth or diffusion-controlled growth 
which will be the aim of further work. Finally, from the 
experimental perspective, it has to be mentioned that the 
model should be able to explain and predict some results. 
Predictions can be made by means of n', co ntroll ed by 
changing the annealing temperature (see Sec. HIE ), and 
c by ion implantation of nucleation centers, or by some 
induced impurities or defects on the sample substrate. 
Some ordered distributions of defects can be induced by 
ion implantation with an appropriate mask, which can 
be trivially introduced in our model. These ideas call for 
further experimental work in order to confirm the validity 
of our model. 
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FIG. 1. Individual grains grown on a square lattice for dif- 
ferent growth probabilities g. 
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FIG. 2. Individual grains grown on a triangular lattice for 
different growth probabilities g. 



a) b) 




FIG. 3. Boundaries between two individual grains obtained 
from simulations with g — 0.5 on a triangular lattice: (a) both 
grains nucleate at the same time, and (b) they nucleate at dif- 
ferent times, yielding a curved interface. 




FIG. 4. log — log plot of the characteristic time tx/2 vs. c in 
the fast nucleation limit over a square lattice: (o) Simulation; 
solid line: power-law fit with slope —0.50 ± 0.01. 




FIG. 5. log — log plot of the characteristic time ti/ 2 '- (a) 
(o) Simulation value, solid line is a power-law fits with slope 
—0.32 ±0.01; (b) Symbols stand for simulation, solid lines are 
power fittings: (o) c = 0.005, slope: 0.34 ± 0.02;(D) c = 0.01, 
slope: 0.34 ± 0.01; (O) c = 0.05, slope: 0.32 ± 0.01; and (A) 
c = 0.1, slope: 0.31 ± 0.02. 
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FIG. 6. Transient KJMA exponent vs. log(t). Circles: 
n — 1 and c = 0.001; squares: n — 0.5 and c = 0.005; di- 
amonds: n = 0.1 and c = 0.05 and triangles: n = 0.01 and 
c = 0.1. g = 1 in all cases. 
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FIG. 7. KJMA exponent dependence on the concentration 
probability c on a 1000 x 1000 triangular lattice. From top to 
bottom: n = 0.001, n = 0.01, n = 0.03, n = 0.07 and n = 1. 
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FIG. 8. KJMA exponent dependence on the nucleation 
probability n on a 1000 x 1000 triangular lattice: (o) 
c = 0.001; (□) c = 0.005; (O) c = 0.01; (A) c = 0.05; and (<) 
c = 0.1, 




FIG. 9. Grain reduced area distribution: (o) Simulation 
with n = 1, g = 1 and c = 0.001; (□) simulation with 
n = 0.001, g = 1 and c = 0.5. Solid line: exact value from 
Eq. (|l|); dashed line: from Eq. @. 
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FIG. 10. Grain reduced area distribution. Simulation with: 
a) n = 0.01 and c = 0.1; b) n = 0.1 and c = 0.05; c) n = 0.5 
and c = 0.005 and d) n — 1 and c = 0.001. g = 1 in all cases. 
Horizontal axis ranges from to 4 and vertical axes from to 
1 in four graphs. 
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FIG. 11. Mean area vs inverse of the potential site concen- FIG. 14. Numerical comparison between normalized distri- 

tration. (o) Simulation values. Dashed line is a linear fit. butions of reduced grain diameter (o), d! , and reduced effec- 

tive diameter, (A') 1/2 (□). 
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FIG. 12. log — log plot of the mean area A: (a) (o) Simu- 
lation value; solid line: power-law fit with slope —0.66 ±0.01; 
(b) (□) simulation, dashed line: power-law fit with slope 
-0.67 ±0.02. 
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FIG. 13. Collapse of the grain effective diameter normal- 
ized distributions for eight g values ranging from 0.05 to 1. 





FIG. 15. Nearest neighbor number normalized distribu- 
tion, (o) n = 1, c = 0.001 (site saturation); (□) n = 0.01 
and c = 0.1 (homogeneous nucleation). Inset: Mean number 
of neighbors, N nn , and its variance, a nn with: a) n = 0.01 
and c = 0.1; b) n = 0.1 and c = 0.05; c) n = 0.5 and c = 0.005 
and d) n = 1 and c = 0.001. 
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FIG. 16. Selfconsistency of Eq. (|l£j) with c = 1 (homo- 
geneous nucleation). (o) Simulation; 
solid line: exponential fit which gives an activation energy 
E a = 2(E n - Eg)/?, = 1.26 ± 0.01, consistent with E„ = 5.1 
eV and E a = 3.2 eV. 







